A new approach to axial coupling constants 

in the QCD sum rule 



Tetsuo NISHIKAWA f]and Sakae SAITO Q 
Department of Physics, Nagoya University, Nagoya 464-8602, Japan 

Yoshihiko KONDOS 

Kokugakuin University, Tokyo 150-8440, Japan 

Abstract 

We derive new QCD sum rules for the axial coupling constants by considering two-point 
correlation functions of the axial-vector currents in a one-nucleon state. The QCD sum 
rules tell us that the axial coupling constants are expressed by nucleon matrix elements of 
quark and gluon operators which are related to the sigma-terms and the moments of parton 
distributions. The results for the iso-vector axial coupling constants and the 8th component 
of the SU(3) octet are in good agreement with experiment. 
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The axial coupling constants are defined by the nucleon matrix elements of the axial- 
vector currents at zero momentum transfer. The iso- vector axial coupling constant Qa 
and the 8th component of the SU(3) octet g^^\o) are known to be expressed by the SU(3) 
parameters F and D: g^^\o) = F + D, gfiO) = 3F - D. In the naive parton model 
the singlet axial coupling constant (/^^(O) is expressed by the fraction of the nucleon spin 
carried by the q quark, Aq, as 5'a'*(0) = + + As, while 5'^''(0) = Am — Ad, and 
9a\^) — + ^d — 2 As. These coupling constants yield an important information on the 
spin structure of the valence and the sea quarks in the nucleon. The parameters F and D are 
precisely known from the measurements in neutron and hyperon (3-decaj experiments, while 
an unexpected small value of (7^'*(0) was found from the EMC data: The quarks contribute 
only a small fraction to the proton's spin, and QaI^), therefore, attracted much attention 

The investigations of the axial coupling constants by QCD sum rules have been done so 
far by the authors in Refs.0-0. Belyaev and Kogan calculated gA (0) by considering a 
two-point correlation function of nucleon currents in an external axial-vector field. loffe 
also calculated (y'^'*(0) and 5'^''(0) with the same correlation function. Their method has some 
difficulties: In the operator product expansion (OPE), there appear new vacuum expectation 
values of quark-gluon composite operators induced by the external field. In addition, the 
spectral function of the correlation function has double and single-pole terms. The residue 
of the former is proportional to (7^(0), and the latter corresponds to the transitions of a 
nucleon state to excited states through the interaction with the external field, but their 
residues are not known. The authors used a method to obtain the values of gA{0) and the 
residue of the single pole term from a fitting procedure, and obtained a good agreement 
with experiment. 

In this paper, we propose a new method to construct QCD sum rules for the axial coupling 
constants from two-point correlation functions of the axial-vector currents in a one-nucleon 
state. With the method, the correlation functions are expressed by nucleon matrix elements 
of quark-gluon composite operators. The spectral function has only a single pole term, 
whose residue is related to (7^(0). As will be shown later, the axial coupling constants 5'^''(0) 
and 5'^'*(0) are related to experimentally or theoretically well-known quantities such as the 
TT-N and K-N sigma-terms and the moments of parton distributions. For the singlet axial 
coupling constant 5'^^(0) we need to fully take into account the chiral anomaly, but it seems 
to be difficult within the ordinary framework of QCD sum rules, and we leave it as a future 
work. 

We first consider a correlation function of axial-vector currents: 

n«(g;P) = ^ J c^xV'?^(T[jg(x),j«(o)])^, (i) 

where the superscript i is the SU{3)f index and = (^jQ)- In Eq-®) the nucleon matrix 
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element is defined by 



(...) 



N 



Y: [{NiPS)\ . . . \N{PS)) - (. . .)o{NiPS)\N{PS))] 



(2) 



where P^ = {E,P) is the nucleon momentum (P^ = M^, M is the nucleon mass), S 
the nucleon spin, (. . .)o the vacuum expectation value, and the one-nucleon state to be 
normalized as {N{PS)\N{P'S')) = {2'Kf6'^{P-P')6ss'- The axial- vector currents are defined 
as 

1 



2^3 



(3) 
(4) 



where u, d and s are the up, down and strange quark fields, respectively. In Eq.(|^), rj^^, = 
Qfj.Qu/ — Qnv is introduced to make the current conserved and suppress, simultaneously, the 
pion contribution to the current 0. 

The Lehmann representation of the correlation function is given by 

,p^^{uj\<^, P) 



n^^(u;,q; P) 



where p^j, is the spectral function defined by 

P/,i.(cj,q;P) = 



duo' 



ijj — ijj' 



TT 



-Imn^,,(u; + ie,q; P). 



(5) 



(6) 



Following Refs.gJ^, we split Xip^y into even and odd parts as Tlp^yiuS) = n^j,(a;^)even + 
ujIlf^y{uj'^)odd, and make a Borel transform on both sides of each Lehmann representation: 



oo 
oo 



even I ~ 

P[n^,(cu2,q; P)odd] = 
where s is a squared Borel mass and B the Borel transformation: 



I duj'uj'exp(-uj''^/s) Pi^,y{uj',ci;P), 
dJ exp {-uj'^l ppy {uj\ q; P) , 



(7) 



B 



lim 

9 

-lij^ — *oo 
n — >oo 
-,,,2 



2\n+l 



d 



(9) 



In Eqs.(l^) and (§) the left hand sides are evaluated by OPE's, which give rise to the Borel 
transformed QCD sum rules. 

Let us now consider the physical content of the spectral function. Among the intermediate 
states of the spectral function, the lowest one is a one-nucleon state. The continuum state 
consists of meson-nucleon states, excited nucleon states and so on. There is an energy gap 
between the pole of the one-nucleon state and the threshold of the continuum states. 
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The contribution of the one-nucleon state to the spectral function is expressed as 



^6{uj + E- V(iV(P^)|j«(0)|iV(P+5'))(iV(P+^')bS(0)|iV(P^)) 



M 



S,S' 



+^6{co -E + E_)l E(^(^^)U2(0)l^(^-^'))(iV(^-^')bg(0)|iV(P^)), (10) 
^- ^ S,S' 

where P± = {E±,P±) = (i? ± cu, P ± q) are the four-momenta of the one-nucleon states. 
From the definitions of the axial coupling constants, the matrix elements in the right hand 
side of Eq. (|T^) are written as 

{NiPS)\jSmNiP'S')) = V,MPS)iXV2)[gfiq')Yl,]uiP'S'), (11) 
(iV(P5)|jg)(0)|iV(P'5')) = uiPS)i\y2)[g^^\q')^,^, + hfiq')q,^,]uiP'S'), (12) 

where u{PS) is a Dirac spinor, A the usual Gell-Mann matrix and q = P' — P. Note that 
the pseudo-coupling term disappears owing to rj^^, in Eq.(|ll]). 

The continuum contribution becomes small in the Borel transformed QCD sum rules, 
since it is exponentially suppressed compared to the one-nucleon state because of the energy- 
gap. Therefore, it is allowed to use a rough model of the continuum: The form of the 
continuum is approximated by the step function with the coefficient being the imaginary 



part of the asymptotic form of the correlation function in the OPE [1C,11]. In the present 
case, however, the continuum contribution to the spectral function is absent within the 
approximation, because from the definition of the correlation function the perturbative part 
is subtracted. This means that the continuum contribution may be very small at least in 
the high energy region. In the following we therefore neglect the continuum contribution to 
the spectral function. 

Hereafter we consider the currents in which the Lorentz indices are contracted, and the 
correlation function in the rest frame of the initial and final nucleon states, P = 0. Thus we 
simplify our notation as follows: n(ci;, q) = U^^^uj, q; M, 0), p{uj, q) = p^^iuj, q; M, 0). Then 
the spectral functions become 



p(3)(o;,q) 



2/ mJm^ + |q|2 



gf{q')\\q'/2-2M' 



X 



5{u + M - + |q|2) _ 5(cj - M + 

2 



M2 + |q|2) 



(13) 



p(«)(a.,q) 



1 



1 



.2^3/ M./M2 + |q|2 



X 



X 



\gf{q^)\\q'/2 - 3M^) + gf{q')hf{q') ■ Mq' - |/.?^(g^)P(gV4) 



5{lu + M- JM^ + |q|2) -S(^-M+ JM^ + |q|2 



(14) 
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Because of crossing symmetry, Eqs.(|TBp and (|T^ are even functions of uj^ so that the Lehmann 
representations of Eq.® is not necessary. From Eqs.(0), (|13D and (0) we obtain 



B 



n(3)| 



1 

M 

,(3)/ 2\|2/ 2 



M2 + |q|2 



exp 



2 



(15) 



1 



1 

M 



1 



M2 + |q| 



M2 + q|2-M Is 



X 



exp 



2.(8)^.2Nr (8)/ 2^ _ (8)/ 2m2 



.(16) 



We expand the right hand sides of Eqs. ([T5|) and (0) in powers of |qp. There is no constant 
term in the contribution of the one-nucleon state to the spectral function. The coefficients 



of |qp in Eqs.(|T3p and (|T^) are proportional to |5'^''(0)p and |(7^''(0)p, respectively. Note 



that /iV(0) contributes to higher order terms because h/A has no singularity at = 0. 
Taking the first derivative with respect to |qp we obtain the desired QCD sum rules at 
12-0: 



|q| 



d 



a|q|2 
d 



5[n(3)(^2^q) 

^(«)(^^q) 



B 



|q|2=0 



|q|2=0 



M 



(17) 



2V3 M 



(8) 



Let us now turn to the OPE of H^^) and n^*^). In the OPE of Eq.([l|), operators of the 
leading terms are of dimension 4. In this work we take into account the terms up to dimension 
6. The result for the iso- vector channel correlation function is in the following: 

■1x2 



n(3)(g) 



N 



2 

^Stto^ (5 (m7mA"w - d^^X'^d) (/x u) 



N 



N 



N 



q=u,d,s 



N 



q6 



q=u,d,! 



N 



^ (i {uS-i^DyDxD„u) ^ + i (dS-i^D^DxD^d 



N 



(19) 
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Similarly, the result for the 8th component of the SU(3) octet is given by 

/ 1 \ 2 



n(«)(g) 



.2^3/ 
Qua. 



10 
7 



{niu (uu) N + nT'd {dd) n + {ss)n 



2g2 \ TT 



a 



N 



Snasq^q'^ 



N 



q" 



N 



3g4 
2'Kasq^q^ 



M7^A"M + J7^A"d + 4s7^A"s) 97^^"?) 

q=u,d,s I 



q=u,d,: 



N 



^ 32q qq q (^uS'y^D^DxD^u) + i ldS-if,D^DxD„d 



N 



N 



(20) 



where -D^'s are covariant derivatives, = GJJ^G"'^'^, and S denotes a symbol which makes 



the operators symmetric and traceless with respect to the Lorentz indices. 

We now discuss about the nucleon matrix elements in Eqs.(0) and (pOl). It is known 
well that mg{qq)i^ is related to the vr-A^ or K-N sigma-term as 



(m„ + md){{uu)N + {dd)N) = 2S^Ar, 
{iris + mu){{ss)N + {uu)n) = '^'^kn- 



(21) 
(22) 



{{as/TT)G'^) AT is expressed by the nucleon mass and mg{qq) jy through the QCD trace anomaly: 



a 
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N 



(23) 



q=u,d,s 

The matrix elements which contain covariant derivatives are related to the parton distribu- 
tions as 



{Sq'y^^D^^---D^„q{fi'^)) 



N 



(24) 



where A„(/i^) is the n-th moment of the parton distributions at scale /x^, and T^^...^^ = 
S [P^-^ ■ ■ ■ Pfin]- For the matrix elements of four quark operators, we apply the factoriza- 
tion hypothesis: In the vacuum, four quark condensates are factorized by the hypoth- 
esis which assumes that the vacuum contribution dominates in the intermediate states: 



(CiC2)o ~ (C^i)o(C^2)o [0)O1- Similarly, for the nucleon matrix elements, we assume 
that the contribution from one nucleon state dominates in the intermediate states P,Q: 
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{Oi02)n ~ {Oi)n{02)q + {Oi)q{02)n- We apply this hypothesis to the following type of 
the nucleon matrix elements, which appear in Eqs.(|l9|) and (pOl): {qf'~i^j.\"'qfqf''~iv\"'qf')N = 
— {S/9)gf^^{qfqf)o{qfqf)]\f6fj', where / and /' are flavor indices. 

We substitute Eqs.(|l9D and into the left hand sides of Eqs. ([T7|) and (p^Sf), respectively. 
Averaging over the iso-spin states, we obtain the QCD sum rules for |5'^''(0)p and |5'^'*(0)p 
as follows: 



9fm 



9fm 



M 
M 



-'ttN 



50 4 

9 9mu + md 



9 ms + rriu 



- + 7{A-{f^') + Ai{f^')) 



A'Kas{qq)i 



352 S 



+ - 



M 

T 

M 



-'■kN 



27 TUu + md 
26 _ 116 m 

y 



+ 



15M3 



3 m„ + rud 



+ 



232 



(25) 



■'KN 



3 nis + mu 



1 iAl{^^') + Ai{^^') + AAf,{^^' 



4'Kas{qq)c 



15M3 r 



352 



■N 



27 niu + nid 



47ras(ss)c 



704 

yr 



2E 



KN 



-'ttN 



rUs +mu mu + rud 



+ 



Alifi^) + Aiifi^) + AAlifi' 



(26) 



respectively, where {qq)o = {uu)q = {dd)^. In Eqs.(p5|) and ( p6D we assume m„ = rUd- From 
these equations, we find that the axial coupling constants are related to the sigma-terms 
and the moments of parton distributions. Since the sigma-terms and the moments are well 



known, we can estimate |fi'^''(0)| and 

We show in Fig. 1 the squared Borel mass s dependence of \gf\'^) \ in Eq.(^ and ^0)1 
in Eq . (p6|) . In plotting the curve in Fig.l, we used the following values of the constants in 
the OPE. The vr- sigma-term is taken from Ref . |]l3l , which are Sttat = 45 MeV. The quark 
masses are taken to be m„ = uid = 7 MeV, = 110 MeV 0. 

Using the above values and the ratio 2{ss) jy / {{uu) jy + {dd)^) = 0.2 given in Ref.[|3 



,(8) 



we can calculate the K-N sigma term averaged over the iso-spin states and the result is 
^KN = 226 MeV. Following Ref. fl^ and adopting the LO scheme in Ref.[0, we calcu- 
lated the moments of parton distributions: A^{1 GeV^) + A^{1 GeV^) = 1.1, A^(l GeV^) + 
A^(lGeV^) = 0.13, ^^(IGeV^) = 0.03, A^(lGeV^) = 0.002. The vacuum condensates are 
taken from Ref.|0|, which are {qq)o = (— 225MeV)^ and {ss)o = 0.8(gg)o. Among the above 
constants in the OPE, the most dominant contribution comes from Al^i^l GeV^)+y42(l GeV^). 

From Fig.l we see good stability of the s dependence for both 5'^^(0) and g^^\o). We see 
that the variations of the curves in a larger s region are small in spite of not including the 
continuum contribution. This fact implies that the continuum contribution is suppressed, 
and the spectral function is allowed to be approximated by the lowest state. 
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s (GeV') 



Fig.l. The squared Borel mass s dependence of i'^ Eq.(|25|) and i'^ Eq. 

The sohd hne corresponds to |5'^''(0)| and the dashed hne to l^'^'*' 



Our estimations of |5'^'*(0)| and |5'^''(0)| taken from the stabihzed region are 

|^??)(0)|^1.2, 1^7^(0)1^0.6. (27) 

The obtained value of |5'a'*(0)| gives a good agreement with the world average (7^''(0) = 
1.260 ±0.002 [|1|]. The value of \gf{0)\ is also close to gf{0) = 0.59 ±0.02 ^ found from 
the data on the baryon octet /5-decays under the assumption of the SU(3) flavor symmetry. 

In summary, we have considered two-point correlation functions of axial-vector currents 
in one-nucleon state and found that the lowest state in the spectral function of the correlation 
function is expressed by the axial coupling constant. We have calculated the iso- vector and 
the 8th component of the SU(3) octet axial coupling constants in the framework of QCD sum 
rules. The results show that the axial coupling constants are expressed by the nucleon matrix 
elements of quark and gluon operators which are related to the sigma-terms and the moments 
of parton distributions. Since the nucleon matrix elements are known well experimentally 
or theoretically, the axial coupling constants are calculated with a small ambiguity, and the 
obtained results are in a good agreement with experiment. 
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Finally, we mention about the singlet axial coupling constant, g^^\o), which is considered 
to be the nucleon spin carried by quarks. The unexpected small value of g^^\o) found by 
EMC has raised a number of understanding of the dynamics of the nucleon spin Within 
the same framework as those for 5'^''(0) and 5'^'*(0) , we find that the calculated ^'^''(0) is 
about 0.8, which is not so small as the value found by EMC. In this calculation, however, 
the effect of the chiral anomaly is not taken into account. 

The authors in Ref. have calculated 5'^^(0) fully taking into account the anomaly 
relation. They have considered a three-point function of nucleon interpolating fields and the 
divergence of the singlet axial- vector current. The form factors, g^^\q'^), are related to the 
vacuum condensates of the quark-gluon composite operators through the double dispersion 
relation. To know at = one must evaluate the correlation function at the zero 



momentum. Although the method to evaluate it is known ||T8[, it involves large uncertainty. 
The calculation of ^'^''(0) in our approach by taking into account the chiral anomaly will be 
reported elsewhere ||19|| . 

The authors (T. N. and Y. K.) would like to thank M. Oka, A. Hosaka and M. Takizawa 
for fruitful discussions. They also wish to acknowledge helpful discussions with O. Morimatsu 
and H. Terazawa. Y. K. thanks the theory group of IPNS(Tanashi) for the hospitality. 
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